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Abstract 



semiconductors, are discussed in the end. 



' We study the pressureless Navier-Stokes-Poisson equations of describing the evolution of 

the gaseous star in astrophysics. The isothermal blowup solutions of Yuen, to the Euler- 
' Poisson equations in , can be extended to the pressureless Navier-Stokes-Poisson equations 

with density-dependent viscosity in _R^. Besides some remarks, about the meaning of the 
QQ ' blowup solutions and the applicability of such solutions to the the drift-diffusion model in 

o ■ 
> 

• 1— I 

^. 

. 1 Introduction 

The evolution of a self-gravitating fluid such as gaseous stars can be formulated by the Navier- 
Stokes-Poisson equation of the following form: 

Pt+V • (pu) =0, 

(pu)j-f V • {pu (g)u) + VP =-pV$ -I- vis{p, u), (1) 
A$(i,a;) = a{N)p, 

where a (TV) is a constant related to the unit ball in i?^: a(l) — 2; a{2) = 2tt and For iV > 3, 



a{N) = N{N - 2)V{N) ^ N{N - 2) 



r(iV/2-fl)' 
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where V{N) is the volume of the unit baU in and F is a Gamma function. And as usual, 
p — p{t,x) and u = u{t,x) 6 are the density, the velocity respectively. P = P{p) is the 
pressure. 

In the above system, the self-gravitational potential field <i> — <i>(t, x) is determined by the 
density p through the Poisson equation. 
And vis(p,u) is the viscosity function: 

vis(p,u) = \7(p{p) V •")■ 

Here we under a common assumption for: 

fi{p) = Kp" 

and K and 9 > are the constants. In particular, when = 0, it returns the expression for the u 
dependent only viscosity function: 

vis{p, u) = kAu. 

The equations and ^2 {vis{p,u) 0) are the compressible Navier-Stokes equations with 
forcing term. The equation ([IJ3 is the Poisson equation through which the gravitational potential 
is determined by the density distribution of the density itself. Thus, we call the system ^ the 
Navier-Stokes-Poisson equations. 

Here, if the vis{p,u) — 0, the system is called the Euler-Poisson equations. In this case, the 
equations can be viewed as a prefect gas model. For = 3, ([T]) is a classical (nonrelativistic) 
description of a galaxy, in astrophysics. See [2], [I] for a detail about the system. 

P = P{p) is the pressure. The 7-law can be applied on the pressure i.e. 

P{p)^Kp^^P^, (2) 
7 

which is a commonly the hypothesis. The constant ^ = cp/c^ > 1, where cp, c„ are the specific 
heats per unit mass under constant pressure and constant volume respectively, is the ratio of the 
specific heats, that is, the adiabatic exponent in ([2]). In particular, the fluid is called isothermal if 
7 = 1. With AT = 0, we call the system is pressureless. 

For the 3-dimensional case, we are interested in the hydrostatic equilibrium specified by u = 
0,5 = In AT. According to 2J, the ratio between the core density p(0) and the mean density p for 
6/5 < 7 < 2 is given by 

P ■ / \ 

-y[z) 



p{0) V ^ 



z=zo 



where y is the solution of the Lane-Emden equation with n — 1/(7— 1), 

y{z) + -y{z) + y[zT = 0, y(0) = a > 0, y(0) = 0, n ^ 



z 7^1 
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and zq is the first zero of y(zo) = 0. We can solve the Lane-Emden equation analytically for 



Vanal{z) = < 



l-^z^, n = 0; 
sin 2; 

, n = 1: 



z ^ 

= , 71 = 5, 



and for the other values, only numerical values can be obtained. It can be shown that for n < 5, 
the radius of polytropic models is finite; for n > 5, the radius is infinite. 

Gambin [4] and Bezard [Ij obtained the existence results about the explicitly stationary solution 
{u = 0) for 7 = 6/5 in Euler-Poisson equations: 

p-(^)''^l+AV)-^/^ (3) 

where A is constant. 

The Poisson equation ([ija can be solved as 



= / G{x-y)pit,y)dy, 

where G is the Green's function for the Poisson equation in the A'^-dimensional spaces defined by 

tl iV = l; 
G{x) = { \og\x\, iV = 2; 

In the following, we always seek solutions in spherical symmetry. Thus, the Poisson equation (fTj)?; 
is transformed to 

r^-^^rr {t, x) + iN- 1) r^-Hr^a (N) pr^-\ 
'"' Jo 

Definition 1 (Blowup) We say a solution blows up if one of the following conditions is satisfied: 

(1) The solution becomes infinitely large at some point x and some finite time Tq; 

(2) The derivative of the solution becomes infinitely large at some point x and some finite time Tq. 

In this paper, we concern about blowup solutions for the 3-dimensional pressureless Navier- 
Stokes-Poisson equations with the density-dependent viscosity, which may describe the phenomenon 
called the core collapsing in evolution of gas star. And our aim is to construct a family of such 
blowup solutions to it. 

Historically in astrophysics, Goldreich and Weber jS] constructed the analytical blowup solution 
(collapsing) of the 3-dimensional Euler-Poisson equation for 7 = 4/3 for the non-rotating gas 
spheres. After that, Makino [7] obtained the rigorously mathematical proof of the existence of 



4 



M.W.Yuen 



such kind of blowup solutions. And in we find tlie extension of the above blowup solutions to 
the case . In [8], the solutions with a from is rewritten as 

For TV > 3 and 7 = (2iV - 2)/iV, 

1 



pit, r) = < 



]vy(^)'^^^'^"'^ forr<a(i)Z^; 



0, for a{t)Z^, < r. 



a(i) 
a{t) 



(4) 



y(z) + 



:y(z)^/(^'2)= fi, yiO) - a > 0, y(0) = 0, 



z - ■ ' {2N -2)K- 
where = [Af(Ar - 2)A]/(27V - 2)K and the finite is the first zero of y{z); 
For TV = 2 and 7=1, 



A 



a{t) 



a{t)= -— ■, a(0) ftp > 0, a(0)= ai; 

a[t} 

1 a(iV) ^("^ 



(5) 



y{x)+-y{x)+ 

X K 



2/(0) = a, 2)(0) = 0, 



where if > 0, = 2\/ K with a sufhciently small A and a are constants. 

For the construction of special analytical solutions to Navicr-Stokes equations, readers may 
refer Yuen's recent results in [TU] . 

In this article, We extend the isothermal B? blowup solutions to the Euler-Poisson equations to 
the pressureless Navier-Stokes-Poisson equations with density-dependent viscosity in i?'^with 9 = \ 
in spherical symmetry: 



Pt + Upr + pUr + -pu = 0, 

r 



p [Ut + UUr) 

in the form of the following theorem. 



47r/9 



pit, s)s^ds + [nphUr + (Kp)(Urr H U, 

r r 



2 , (6) 



Theorem 2 For the 3 -dimensional pressureless Navier-Stokes-Poisson equations with 9 — 1, in 
spherical symmetry, there exists a family of solutions, 

1 .r.urr . . -C 

(7) 



Pit,r) = 



T-Ct ' 



X kC 

where T > 0, k > 0, C>0 and a are constants. 



(T - Ci)3 

2 An ^(^^ 

y{x) + -y{xH—e = 0, y{0) = a, y{0) = 0, 



And the solutions blowup in the finite time T/C . 



2 Separable Blowup Solutions 

In this section, before presenting the proof of Theorem [2l we prepare some lemmas. First, we 
obtain a general class of solutions for the continuity equation of mass in radial symmetry ©i. 
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Lemma 3 For the 3- dimensional conservation of mass in spherical symmetry 

2 

pt + upr + pUr + -pu = 0, (8) 
r 

there exist solutions, 

where T and C are positive constants and f £ is a non-negative function . 

Proof. We just plug © into dH]). Then 
2 

Pt + Upr + pUr H pu 

r 

_ (_3)(_C)e/('-/(r-g*)) e-^('-/(^-c*))/(r/(r - Ct)) r{-l){-C) 
{T-Ct)'^ [T-Ct)^ [T-Cty 

i-C) e/('^/(^-^*))/(r/(r-Ci)) 1 e-^f'-Zt^-*)) {-C) 2 e-^('-/(^-C't)) (_^) 



T-Ct (T - Ct)3 T - Ct (T - Ct)-'' T-Ct r (T - Ct)3 T - Ci 

_ 3Ce-^('-/(^-^*)) Ce/('-/(^-^*))/(r/(T - Ct)) 
{T-CtY ^ {T-Ctf ^ 

Ce/('-/(^-g«))/(r/(r-Ct)) Ce^('^/(^-g«)) 2Ce^('-/(^-g*)) 
(T - Ct)5 (T - Ct)4 (T - Ct)4 

= 0. 

The proof is completed. ■ 

Besides, we need the two lemmas for stating the property of the function y{z). The similar 
lemmas are already given in (Lemmas 9 and 10) Yuen's article [9]. For the completeness of the 
article, the proof of the below lemmas are presented. 

Lemma 4 There exists a sufficiently small xq > 0, such that the equation 

j/(a:)+^y(a:)+ae^(-)=0, ^^^^ 
2/(0) = a, 2/(0) = 0, 

where a > and a are constants, has a solution y — y{x) e C^[0, xq\. 

Proof. The lemma can be proved by the fixed point theorem. Multiply (fTU)) by x, to give 



d 
dx 

Notice y(0) = 0, we have 



{x^y{x)) 



9 
X^ 



y{x) = -J- / aey'^'h'^ds. 
By using 2/(0) — a, pH)) is reduced to 

y{x) = ^ / ae^f^^s^ds, 1/(0) = a. 
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Set 



Jo 

then for any xq > 0, we get / € C^[0, x^]. and for any yi^ 7/2 & C^[0, xo], we have, 



\f{x,yi{x)) - f{x,y2{x))\ = 



As is a C"'^ function of y, we can show that the function e^, is Lipschitz-continuous. And we get, 

\f{x,yi{x)) - f{x,y2{x)\ 

_ 0{l)f^s^\{y2{s)-yi{s)\ds 



where r G [0, a;] C [0, a;o]. Let 



< O(l)a;o sup |?/i(s) - y2(s)l 

0<X<Xn 



Ty{x)^a+I f{s,y{s))ds. 
Jo 



We have Ty € C[0, xq] and 



\Tyi{x) - Ty2ix)\ 



f{s,yi{s))ds - f{s,y2is))ds 
Jo 



< O(l)a;o sup \y{x)i - y{x)2\ 

0<x<~xo 



By choosing ccq > to be a sufficiently small number, such that 0{1)xq < 1, this shows that the 
mapping T : C[0, Xq] C[0, xq], is a contraction with the sup-norm. By the fixed point theorem, 
there exists a unique y{x) G C[0, xq], such that Ty{x) — y{x). The proof is completed. ■ 



Lemma 5 The ODE, 

2 

j/(a;)+-y(x)+CTe^("')= 0, 

X 

2/(0) = a, 2/(0) = 0, 
where a > ari(i a are constants, has a solution in [0, +00) and lim y{x) 

x^-\-oo 



(11) 



-00. 



Proof. By integrating pT|) . we have. 



(12) 



Thus, for Q < X < Xq, y{x) has a uniform lower upper bound 



y{x) < 2/(0) = a. 



As we obtained he local existence in Lemma IH there are two possibilities: 
(1)2/(2;) only exists in some finite interval [0,a;o]: 



Analytical Blowup Solutions 



7 



(la) lim y{x) — — oo; 

(Ib)y(x) has an uniformly lower bound, i.e. y{x) > ao for some constant ao- 
(2)y{x) exists in [0, +00): 
(2a) lim y{x) — —oo\ 

X — * + oo 

{2\))y{x) has an uniformly lower bound, i.e. y{x) > ao for some constant ao. 

We claim that possibihty (1) doesn't exist. We need to reject (lb) first: If the statement (lb) is 
true, lfT2|) becomes 

= re-^s^dsKyix). (13) 

■J X Jo 

Thus, y{x) is bounded in [0,a;o]. Therefore, we can use the fixed point theorem again to obtain a 
large domain of existence, such that [0, xq + iJ] for some positive number S. There is a contradiction. 
Therefore, (lb) is rejected. 

Next, we don't accept (la) because of the following reason: It is impossible that lim y{x) = —00, 

X — >xo — 

as from (|13p . y{x) has a lower bound in [0, xq]: 

crxQe'^' 



< V{x). (14) 



Thus, (dU becomes. 



rxo 

y{xo) = y(0) + / yix)dx 



axe" , 
> a — — - — ax 







2„a 



a 



ax^e 



Since y{x) is bounded below in [0, a;o], it contracts the statement (la), such that lim y{x) = —00. 

X — >xo- 

So, we can exclude the possibility (1). 

We claim that the possibility (2b) doesn't exist. It is because 



Jo a;^ Jo 



Then, we have. 



ere"" 

y{x) <a —x^. (15) 



By letting x 00, (|T5| turns out to be, 

y(x) = -00. 



Since a contradiction is established, we exclude the possibility (2b). Thus, the equation pTjl exists 
in [0, +00) and lim y{x) = —00. This completes the proof. ■ 

x — ^+00 

Up to this stage, we are already to give the proof of Theorem [2] 
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Proof of Theorem 2. From Lemma [31 we easily get that (O satisfy (jS])!. For the momentum 
equation (02, we get, 



p(ut + UUr) ^ 5- / p(t, s)s^ds - [p{p)]rUr — p{p){Urr + -Ur ^-u) 



(T - cty 
i-c) 

T-Ct 
^Cey^r/iT-ct)) 



T-Ct T -Ct 



Anp 



r2 J {T-Ct)^ 



s ds 



{Kp)r 



Pip) 



2 (-1) 2 (-1) 



T-Ct r^T-Ct 



1 C Anp /•e?^(-/(^-")) 



(r-Ct)3 "^T-Ct' T-CtT-Ct J (T-Ct) 



-s^ds 



kCp r Attp /•e?^(^/(7'-«)) 



[T -Ctf^T -Ct' r2 y (T-Ct) 



nCiji 



An 



T-Ct' r^{T-Ct) 



By letting w = s/(T - Ct), it follows 



(T - Ct)2 



^ • / "f" N 47r 



T/{T-Ct) 



T-Ct (yzcl)^ 



[T-Ctf^ \T-CtJ ' 
And denote z — r/{T — ct), 



Differentiate Q[z) with respect to 2, 

Q(z) - KCy{z)+ATTey^'-^ + ^"^^ " / ey^^^Lo'^ds 



(16) 



(17) 



(18) 



2.CC. , , (-2)47r /• ^(^) 
= y(z) H ^ — / e'''' du) 

z Z'^ J 



= --Q(^), 

z 

where the above result is due to the fact that we choose the following ordinary differential equation: 

2 47r ^(^^ 

y{z) + -y{z)+—e =0 
z kC 

2/(0) = a, m= 0. 

With Q(0) — 0, this implies that Q{z) — 0. Thus, the momentum equation ^2 is satisfied. 
With the Lemma [5] about y{z), we are able to show that the family of the solutions blows up in 
finite time T/C. This completes the proof. ■ 
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The statement about the blowup rate will be immediately followed: 

Corollary 6 The blowup rate of the solution ^ is 

^jim^p(t,0)(T-Ct)3>O(l). 

Given that the sign of the constant C in ([5]) is changed to be negative, the below corollary is 
clearly shown. 

Corollary 7 For the 3 -dimensional pressureless Navier-Stokes-Poisson equations with 6 — 1 in 
spherical symmetry, 0), there exists a family of solutions, 

i'^-^'t (19) 
y{x) + -y{x) +—e = 0, y(0) = a, y{Q) = 0, 
X kC 

where T > 0, k > 0, C < and a are constants. 

At last, we give some remarks to discuss the constructed blowup solutions. 

Remark 8 Our solutions are the modification versions of blowup solutions of Yuen in the 
2-dimensional Euler-Poisson equations. However, we emphasize that the process is not trivial, as 
the blowup solutions cannot be generated to the other dimensional case {N ^ 2). And the blowup 
solutions ^ of Goldreich, Weber and Makino also cannot be generated to the lower dimensions '(1- 
dimensional or 2-dimensional) case. 

Remark 9 It is well-known that it is an open problem to construct the analytical blowup solutions 

in the i- dimensional case to understand the blowup phenomena (black hole) in the Euler-Poisson 

equations in astrophysics. But for the non-trivial blowup solutions, there only exists the blowup 

solutions Meanwhile, due to the solution structure of u in ^ and 

a{t) 
u = —r-r, 
a{t) 

where a{t) ^ is a function of time, 

that automatically satisfies that the density-independent viscosity function {6 — Q): 

vis(u) — kAu — K[Urr H Ur ^ U) = 0. 

In this article, we shift the scope of consideration from the Euler-Poisson equations in to the 
Navier-Stokes-Poisson equations with density- dependent viscosity. We can extend the isothermal 
blowup solutions to the Euler-Poisson equations to pressureless Navier-Stokes-Poisson equations 
with density- dependent viscosity in . Here blowup solutions can answer the above problem in 
the modified i- dimensional fluid dynamic system. We notice that the extension of such blowup 
solutions are not suitable for solutions from R^ to R^^^ (N > 3J. In the future, the blowup 
solutions in the Euler-Poisson equations will still be sought. 
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Remark 10 Besides, if we consider the 'i- dimensional drift- diffusion model in semiconductors in 
spherical symmetry, 

2 

pt + UPr + pUr + -pu = 0, 

r ^ 2 2 

p {Ut + UUr) = Jq p{t, S)s'^ds + [Kp]rUr + {Kp){Urr + -Ur - -^u), 

the special solutions with infinite mass may be obtained as follows: 



nit r) = - e''(''/('^-'^*)) u(t r) = —r- 

' ' {T-Ctf ' ^ ' ' T-Ct' 

y{x)+-y{x) - — e = 0, 2/(0) = a, y(0) = 

X 



(20) 
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